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Abstract: We consider networks of massive particles connected by non-linear springs. 
Some particles interact with heat baths at different temperatures, which are modeled 
as stochastic driving forces. The structure of the network is arbitrary, but the motion 
of each particle is ID. For polynomial interactions, we give sufficient conditions for 
Hormander's "bracket condition" to hold, which implies the uniqueness of the steady 
state (if it exists), as well as the controllability of the associated system in control theory. 
These conditions are constructive; they are formulated in terms of inequivalence of the 
forces (modulo translations) and/or conditions on the topology of the connections. We 
illustrate our results with examples, including "conducting chains" of variable cross- 
section. This then extends the results for a simple chain obtained in [5] . 
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1. Introduction 

We consider a network of interacting particles described by an undirected graph Q = 
(V, £) with a set V of vertices and a set £ of edges. Each vertex represents a particle, 
and each edge represents a spring connecting two particles. We single out a set V* C V 



2 



N. Cuneo and J. -P. Eckmann 



of particles, each of which interacts with a heat bath. We address the question of when 
such a system has a unique stationary state. This question has been studied for several 
special cases: Starting from a linear chain [5,4], results have become more refined in 
terms of the relation between the spring potentials and the pinning potentials which tie 
the masses to the laboratory frame [11,2]. This problem is very delicate, as is apparent 
from the extensive study in [8] for the case of only 2 masses. 

We provide conditions on the interaction potentials that imply Hormander's "bracket 
condition," from which it follows that the semigroup associated to the process has a 
smoothing effect. This, together with some stability assumptions, implies the unique- 
ness of the stationary state. The existence is not discussed in this paper, but seems well 
understood in the sense that the interaction potentials must be somehow stronger than 
the pinning potentials. This issue will be explained in a forthcoming paper [3]. 

Since the problem is known (see for example [7] and [1]) to be tightly related to the 
control problem where the stochastic driving forces are replaced with deterministic con- 
trol forces, we shall use the terminology of control theory, and mention the implications 
of our results from the control-theoretic viewpoint. 

We work with unit masses and interaction potentials that are polynomials of de- 
gree at least 3, and we say that two such potentials V\ and Vi have equivalent second 
derivative if there is a S € R such that V{'( • ) = V^'i ■ +5). 

We start with the set V* of particles that interact with heat baths, and are therefore 
controllable. One of our results (Corollary 5.6) is formulated as a condition for some of 
the particles in the set of first neighbors W(V») of V* to be also controllable. Basically, 
the condition is that these particles must be "inequivalent" in a sense that involves both 
the topology of their connections to V* and the corresponding interaction potentials. 
More precisely, a sufficient condition for a particle v e 7\A(V») to be controllable is that 
for each other particle w e W(V*) at least one of the two conditions holds: 

(a) v and w are connected to V» in a topologically different way, 

(b) there is a particle c in V* such that the interaction potential between c and v and 
that between c and w have inequivalent second derivative. 

It is then possible to use this condition recursively, taking control of more and more 
masses at each step (Theorem 5.7). If by doing so we can control all the masses in the 
graph, then Hormander's bracket condition holds. 

In Sect. 6 we give examples of physically relevant networks whose controllability 
can be established using this method. 

Our results imply in particular that connected graphs are controllable for "almost 
all" choices of the interaction potentials, provided that they are polynomials of degree 
at least 3 (Corollary 6.3). 

2. The system 

We define a Hamiltonian for the graph Q = (V, £ ) as follows. Each particle v e V has 
position q v 6 R and momentum p„eR and is "pinned down" by a potential U v (q v ). 
Throughout, we assume the masses being 1, for simplicity of notation. See Remark 4.20 
on how to adapt the results when the masses are not all equal. 

We denote each edge e 6 £ by {u, v} (or equivalently {v, u}) where u, v are the 
vertices adjacent to e. 1 To each edge e = {u, v}, we associate an interaction potential 

1 Due to the physical nature of the problem, we assume that the graph has no self-edge. 
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V uv (q u - q v ), or equivalently V vu (q v - q u ) with 

V vu (q v - q u ) = V uv (q u - q v ) . (2.1) 

Note that we do not require the potentials to be even functions; the condition (2.1) 
just makes sure that considering e = {u, v} or e = {v, u} leads to the same physical 
interaction, which is consistent with the fact that the edges are not oriented. 
With the notation q = (q v ) ve v an d P = (Pv)vev the Hamiltonian is then 

H(q,p) = ]T (p 2 j2 + U v (q v )) + ]T V e (Sq e ) , 

where it is understood that V e (Sq e ) denotes one of the two expressions in (2.1). 
Finally, we make the following assumptions: 

Assumption 2.1 

1. All functions are smooth. 

2. The level sets of H are compact, i.e., for each K > the set {(q,p) \ H{q 1 p) < K} 
is compact. 

3. The function exp(—/3H) is integrable for some (3 > 0. 

Each particle v G V* is coupled to a heat bath at temperature T v > with coupling 
constant j v > 0. For convenience, we set j v = when v V* . The model is then 
described by the system of stochastic differential equations 

dq v =p v dt , (2.2) 
dp v = -U' v {q v )dt - 8 qv (J^ V e {5q e ))dt - lvPv dt + y / 2T vlv dW v (t) , 

where the W v are identical independent Wiener processes. The solutions to (2.2) form 
a Markov process. The generator of the associated semigroup is given by 

L = X ° + E ^T v d 2 pv , 
vev, 

with 

= - 1vPvd Pv + Y (Pvdqv - K(lv)d Pv ) - Y Kvilu - q v ) ■ (d Pu - d Pv ) . 

vGV, veV {u,v}e£ 

From now on, we assume that the interaction potentials V e , e e £ are polynomials of 
degree at least 3. The condition on the degree means that we require throughout the 
presence of non-harmonicities. The fully-harmonic case has been described earlier [6], 
and the case where some but not all the potentials are harmonic is not covered here. 
We will show in a counter-example (Example 6.8) that the non-harmonicities are really 
essential for our results. We make no assumption about the pinning potentials U v \ we 
do not require them to be polynomials, and some or all of them may be identically zero. 

We work in the space R 2 I V I with coordinates x = (q,p). We identify the vector 
fields over R 2 I V I and the corresponding first-order differential operators in the usual 
way. This enables us to consider Lie algebras of vector fields over R 2 ' v l, where the Lie 
bracket [ • , • ] is the usual commutator of two operators. 
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Definition 2.2. We define M. as the smallest Lie algebra that 

(i) contains d Pv for all v G V*, 2 

(ii) is closed under the operation [ ■ , Xq], 

( Hi) is closed under multiplication by smooth scalar functions. 

By the definition of a Lie algebra, M is closed under linear combinations and Lie 
brackets. 

Definition 2.3. We say that a particle v G V is controllable if we have d qv , d Pv G A4. 
We say that the network Q is controllable if all the particles are controllable, i.e., if 

3,„a p ,eM forallveV. (2.3) 

The aim of this paper is to give sufficient conditions on Q and the interaction poten- 
tials, which guarantee that the network is controllable. 

If the network is controllable in the sense (2.3), then Hormander's condition 3 [9] 
holds: for all x, the vector fields F G M evaluated at x span all of R 2 I V I, i.e., 

{F(x) | F G M} = R 2|v| for all x G R 2|v| . (2.4) 

Hormander's condition implies that the transition probabilities of the Markov pro- 
cess (2.2) are smooth, and that so is any invariant measure (see for example [1, Cor. 7.2]). 
We now briefly mention two implications of these smoothness properties. Proposi- 
tion 2.4 and Proposition 2.5 below can be deduced from arguments similar to those 
exposed in [7], and will be discussed in more detail in the forthcoming paper [3]. 

Proposition 2.4. Under Assumption 2.1, if (2.4) holds, then the Markov process (2.2) 
has at most one invariant probability measure. 

The control-theoretic problem corresponding to (2.2) is the system of ordinary dif- 
ferential equations 

q v =p v , (2.5) 
Pv = -U' v (q v ) - d qv Ve(£<7e)) + ( u v(t) - JvPv) ■ Uev, , 

e€£ 

where for each v G V*, u v : R — >• R is a smooth control function (i.e., the stochastic 
driving forces have been replaced with deterministic functions). 4 

Proposition 2.5. Under the hypotheses of Proposition 2.4, the system (2.5) is control- 
lable in the sense that for each x^ = (q^ ,p^) and x^ — (q^\p^), there are a 
time T and some smooth controls u v , v G V*, such that the solution x(t) of (2.5) with 
x(0) = verifies x(T) — x^\ 

In fact, (2.4) is a well-known condition in control theory. See for example [10], 
which addresses the case of piecewise constant control functions. In particular, (2.4) 
implies by [10, Thm. 3.3] that for every initial condition a;' ' and each time T > 0, 
the set A(x^ , T) of all points that are accessible at time T (by choosing appropriate 
controls) is connected and full-dimensional. 

2 Due to the identification mentioned above, we view here c9 PtJ as a constant vector field over R. 2 v I . 

3 The condition (2.4) is slightly different, but equivalent to the usual statement of Hormander's criterion. 
This can be checked easily. In particular, closing M under multiplication by smooth scalar functions does 
not alter the set in (2.4), and will be very convenient. 

4 Whether or not we keep the dissipative terms — fvfv in (2.5) makes no difference since they can always 
be absorbed in the control functions. 
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3. Strategy 

We want to show that d Qv , d Pv G M. for all ueV. The next lemma shows that we only 
need to worry about the d Pv . 

Lemma 3.1. Let Abe a subset of V. 

V J2 d P- eM then E d i- e M ■ 

veA veA 

Proof. Assuming J2 v eA ®Pv e M, we ^ n( ^ tnat 

[ e d P, >*o] = E ^ - E ^ (3 -v 

veA veA vev,nA 

is in M. But since d Pv E M for all v e V*, the linear structure of M implies 

J2vev nA 1vd Pv £ M. Adding this to the vector field in (3.1) shows that 

J2veA dqv e ^> as claimed. □ 

Definition 3.2. We say that a set A C V /s jointly controllable i/5Z„ gy i is in M 
(and therefore, also ^2 veA d qv by Lemma 3.1). 

Requiring all the particles in a set A to be (individually) controllable is stronger than 
asking the set A to be jointly controllable (indeed, if all the d Pv , v e A are in M, then 
so is their sum). We will obtain jointly controllable sets and then "refine" them until we 
control particles individually. 

The strategy is as follows. In the next section, we start with a controllable particle c, 
and show that its first neighbors split into jointly controllable sets. Then, in Sect. 5, we 
consider several controllable particles, and basically intersect the jointly controllable 
sets obtained for each of them, in order to control "new" particles individually. Finally, 
we iterate this procedure, taking control of more particles at each step, until we establish 
(under some conditions) the controllability of the whole network. 

Remark 3.3. Observe in the following that our results neither involve the interaction 
potentials U v nor the coupling constants -f v . 

4. The neighbors of one controllable particle 

We consider in this section a particle c e V, and denote by T c the set of its first neigh- 
bors (the "targets"). The following notion of equivalence among polynomials enables 
us to split T c into equivalence classes. 

Definition 4.1. Two polynomials f and g are called equivalent if there isaieR such 
that /(•) =g( - +6). 

Definition 4.2. We say that two particles v,u G T c are equivalent (with respect to c) if 
the two polynomials V c " and V" u are equivalent. 

Since this relation is symmetric and transitive, the set T c is naturally decomposed 
into a disjoint union of equivalence classes: 

r c - u 4 7; c . 

An explanation of why we use the second derivative of the potentials instead of the first 
one (i.e., the force) will be given in Example 6.7. The main result of this section is 
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Theorem 4.3. Assume that c is controllable. Then, each equivalence class 7~ c is jointly 
controllable, i.e., 

d Pv EM for alii. (4.1) 

Furthermore, there are constants S cv such that for all i, 

(ic - q v + 5 cv )d Pv EM. (4.2) 

The second part of the theorem will be used in the next section to intersect the 
equivalence classes of several controllable particles c. We will now prepare the 
proof of Theorem 4.3. We assume in the remainder of this section that c is controllable. 
And since c is fixed, we write T and % instead of T c and T£. 

Lemma 4.4. We have 

J2V^(q c -q v )d Pv EM. (4.3) 

Proof. From Lemma 3.1 we conclude that d qc E M. Therefore, we find that 

[d qe ,X ] = -K'(q c )d Pc £ V£(q c - q v ) ■ (d Pc - d Pv ) 

is in Ai. Now, since d Pc E Ai and since M. is closed under multiplication by scalar 
functions, we can subtract all the contributions that are along d Pc and obtain (4.3). □ 

We need a bit of technology to deal with equivalent polynomials. 

Definition 4.5. Let g(t) — X^j=o ai ^ % A' be a polynomial of degree k > 1. Ifa^-i = 0, 
we say that g is adjusted. As can be checked, the polynomial g(-) = g( - — ak-i/a,k) 
is always adjusted, and is referred to as the adjusted representation of g. 

Observe that a polynomial and its adjusted representation are by construction equiv- 
alent and have the same degree and the same leading coefficient. In fact, given a poly- 
nomial g of degree k > 1, g is the only polynomial equivalent to g that is adjusted. This 
adjusted representation will prove to be very useful thanks to the following obvious 

Lemma 4.6. Two polynomials g and h of degree at least 1 are equivalent iff g = h. 

Remark 4.7. If all the interaction potentials are even, then all the V c " are automatically 
adjusted, and some parts of the following discussion can be simplified. 

We shift the argument of each V" v by a constant S cv so that they are all adjusted. We 
let /„ be the adjusted representation of V c " and use the notation 

<^v — qc qv &cv 

so that 

fv {x v ) = V"v (qc- qv) for all q c , q v e R . 
With this notation, (4.3) reads as 

fv(*v)d Pv G M . (4.4) 

v£T 

We will now mostly deal with "diagonal" vector fields, i.e., vector fields of the kind 
(4.4), where the component along d Pv depends only on x v . When taking commutators, 
it is crucial to remember that x v is only a notation for q c — q v + S cv . 
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Remark 4.8. By the definition of equivalence and Lemma 4.6, two particles v,w G T 
are equivalent iff f v and f w coincide. 

Lemma 4.9. Consider some functions g v , v G T. 

Proof. This is immediate by commuting with d qc (which is in M by Lemma 3.1). □ 
We now introduce the main tool. 

Definition 4.10. Given two vector fields Y and Z, we define the double commutator 

[y : Z] by 

lY:Z}=[[X ,Y\,Z\- 
Obviously, if the vector fields Y and Z are in M., then so is [y : Zj . 
Lemma 4.11. Consider some functions g v , h v , v G T. 7/ien 

[X! 9v(x v )d Pv : h v(xv')d Pv ,] = ^2(g v h v )' \x v ) d Pv . (4.6) 

Proof. Observe first that (omitting the arguments x v ) 

[Xo, ^2gvd Pv ] = ^2(p c - Pv)g' v d Pv - ^ 9v9 qv + ~/v9vd Pv . 
veT vgT vgT veT nv„ 

Commuting with J2v'eT ( Xv ' )®p v > gi yes the desired result. □ 

We will prove Theorem 4.3 starting from (4.4) and using only (4.5) and double 
commutators of the kind (4.6). 

Let d v be the degree of /„. Note that since the interaction potentials are of degree at 
least 3, we have d v > 1. We define 

d = max d v > 1 

veT 

as the maximal degree of the adjusted polynomials /„ with veT. We can then write 

d 

fv{x) = J2 b ^ xJ /j- ! 

3=0 

for some real coefficients b v j, j = 0, . . . ,d, with 

b vj = if j > d v and b v , dv -\ = , 

for all veT. 

Definition 4.12. We define the set of particles veT corresponding to the maximal 
degree d: 

V d = {v eT\d v = d} . 
For every I, < I < d, w e define the set 

Bf = {b vl \veV d } 

of distinct values taken by the coefficients ofx^/t\ in f v , v G T> d . 
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We begin with a technical lemma. Observe how it is expressed in terms of the x v . In a 
sense, this shows that the x v are really the "natural" variables for this problem. Thus, 
in addition to making the notion of equivalence trivial (Remark 4.8), working with 
adjusted representations will be very convenient from a technical point of view. 

Lemma 4.13. The following hold: 

( i) For each b € Bf , we have 

^2 x v d Pv G M and ^ d Pv £ M . (4.7) 

veV d :b vd =b vEV d :b vd =b 

( ii) Furthermore, 

^2 x v9 Pv G M and ^ 9 Pv € M . (4.8) 

vev d vev d 

(Hi) Let a v , f3 v , v <E T) d be real constants. If d > 2, we have the two implications 

if a v d Pv G M , then a v x v d Pv G M , (4.9) 

veT> d vev d 

ifJ2 a " d Pv> J2P vd Pv eM > then ^2 <*v Pv d Pv G M . (4.10) 

veT> d vET> d v£T> d 

Remark 4.14. Observe that the assumption d > 1 is crucial in the proof of (i). Requiring 
the /„ to be non-constant ensures that we can find non-trivial double commutators, 
which is the crux of our analysis. See Example 6.8 for what goes wrong for harmonic 
potentials. 

Proof, (i). By (4.4) and using (4.5) recursively d — 1 times, we find that 

Y = ^2(d d ~ 1 fv)(x v )d Pv = ( b v,d-i +bvdX v )d Pv 

is in M. But now, by (4.6), 

\Y : Y/2] = bvdiKd-i + b vd x v )d Pv G M . 

Taking more double commutators with Y/2, we obtain for all r > 1: 

^2 Kd ( b v,d-i + KdX v )d Pv E M . 

But the sum above is really only over D d since b v d ^ only if v G "D . Moreover, for 
these v, we have b Vi d-i = since the polynomials are adjusted, so that for all i > 2, 

^'„M,^' (4.H) 
vev d 
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Let b G B d . Using Lemma A.l with s = 1 and with the set of distinct and non-zero 
values {b v d | v G T> d } = B d we find real numbers ri , r2 , . . . , r n (with n = | Bf | ) such 
that Y^i=i r i ^ % vd e q ua l s 1 if b vd = b and when b vd ^ b. Thus, 

n 

r « 12 o s »^p» = 5Z x vd Pv 

i=l veV d veV d :Kd=b 

is in M by (4.1 1). This together with (4.5) establishes the second inclusion of (4.7), so 
that we have shown (i). 

The statement (ii) follows by summing (i) over all b G Bf. 

Proof of (iii). Let us assume that d > 2 and that J2vev d a v ®p v £ -A4 . By (4.7), we 
have for each b G Bf that 

\ x ^ = 12 yfv* e M ■ 

v£V d : b vd =b veV d :b vd =b V 

Taking the double commutator with J2vev d a « ®Pv an d summing over all b G B d 
shows that 




Since we assume here d > 2, we have Z = ^2 ver {d d 2 f v ) {x v )d Pll G M. But then, 




is also in M. Recalling that b v . d - X = for all v G V d , we obtain (4.9). Finally (4.10) 
follows from (4.9) and the double commutator 

lY a v x v d Pv : Y Pv d Pv j = ^ a v f3 v d Pv . 

vev d vev d vev d 

This completes the proof. □ 

With these preparations, we can now prove Theorem 4.3. 

Proof (of Theorem 4.3). 

We distinguish the cases d = 1 and d > 2. 
Case d = 1: This case is easy. Since all the /„ have degree 1, we have that f v (x v ) = 
b v \ x v for all v G T, with b vl ^ 0. Consequently, the sets % consist of those v which 
have the same b v \ (see Remark 4.8). Thus, we have by (4.7) that for each 

^3 P ,£^ and x " d v, e ^ (if = 1) . (4.12) 

veTi v£% 

This shows that the conclusion of Theorem 4.3 holds when d = 1. 
Case d > 2: In this case, (4.7) is not enough. First, (4.7) says nothing about the masses 
v G T \ V d , for which 6^ = 0. And second, (4.7) provides us with no way to "split" 
the d Pv corresponding to a common (non-zero) value b of b v d, even though the corre- 
sponding v might be inequivalent due to some b v k with k < d. To fully make use of 
these coefficients, we must develop some more advanced machinery. 
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Definition 4.15. We denote by Vd the vector space of real polynomials in one variable 
of degree at most d. We consider the operator G : Vd — > Vd defined by 

(Gv)(x) = (x-v(x))' , 

and we introduce the set of operators 

T = sp&n{G,G 2 ,...,G d+1 } . 

Observe that by (4.4) and (4.8) we have 

1^2 fv{x v )dp v : ^2 x v d Pv ] = ^2 ( G fv)(x v )d Pv € M . 
veT vev d vev d 

Note that we obtain a sum over V d only. By taking more double commutators with 
J2 V £V d x vd Pv , we find that J2vev d (G k fv)(xv)d Pv is in M for all k > 1. Thus, by the 
linear structure of M, we obtain 

Lemma 4.16. For all P e T, we have 

Y2, (Vfv)(x v )d Pv e M . 

vev d 

It is crucial to understand that it is the same operator P that is applied simultaneously 
to all the components, and that the components in T \ T> d are "projected out." 
We now show that some very useful operators are in T . 

Proposition 4.17. The following hold: 

( i) The projector 

d 

Sr.Vd^Vd, Yl bixi / i] ^ bext / & 

i=0 

belongs to T for all I = 0, . . . ,d. 

( ii) The identity operator 1 acting on Vd is in T. 

Proof. Consider the basis B = (eo, e\, . . . , e^) of Vd where ej(x) = x' J j ' j\ . Observe 
that for all j > we have Gej = (j + l)ej, so that G is diagonal in the basis B. Thus, 
G fe is represented by the matrix diag(l' £ , 2 k , . . . , (d+ l) k ) for all k > 1. Consequently, 
for each t G {0, . . . , d}, there is by Lemma A.l with s = a linear combination of 
G, G 2 , . . . , G d+1 that is equal to S^. This proves (i). Moreover, we have that J2i=o = 
1, so that the proof of (ii) is complete. □ 

Lemma 4.18. For all £ = 0, . . . , d, and for each b E Bf = {b vg \ v e V d } we have 



(4.13) 
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Proof. Let I G {0, 1, . . . , d}. Using Lemma 4.16 and Proposition 4.17(i) we find that 
J2vevd( b vex e /i\)d Pv is in M. Using (4.5) repeatedly, we find that J2vev d b ved Pv is 
in M. Thus, by (4. 10), 

Ke d Pv e M for alii > 1. 

vev d 

Then, applying Lemma A. 1 to the set Bf \ {0} and with ,s = 0, we conclude that 

d Pv eM for all b G Bf \ {0} . (4.14) 

veV d : b vt =b 

If ^ Bf, we are done. Else, we obtain that (4.13) holds also for b = by summing 
the vector field (4. 14) over all b G Bf\ {0} and subtracting the result from J2vev d ®Pv 
(which is in M by (4.8)). This completes the proof. □ 

Remember that by Remark 4.8, a given equivalence class 77 is either a subset of V d 
or completely disjoint from it. 

Lemma 4.19. Let % be an equivalence class such that % C V d . Then 

<9 P „ G M , and ^ Xr " d Pv e • (4.15) 

ve% veTi 

Proof. All the polynomials /„, v G % are equal. Thus, there are coefficients q 6 Bf, 
£ = 0,l,...,d such that 

d 

%= f){veV d \b ve = c e } . (4.16) 

e=o 

By Lemma 4.18, we have for alH = 0, . . . , d that 

]T d Pv G M . (4.17) 

Now observe that whenever two sets B,B' C 2? d are such that ^2 vE b 9 Pv G 
and X)«eB' 9 P „ G .M, we have by (4.10) that J2veBnB> ®Pv <= M. Applying this 
recursively to the intersection in (4.16) and using (4.17) shows that X^er ®Pv e 
Using now (4.9) implies that J2 v eTi x vdp v e -M> which completes the proof. □ 

With these tools, we are now ready to complete the proof of Theorem 4.3 (for the 
case d > 2). By Lemma 4.19, we are done if V d = T (i.e., if all the f v ,veT have 
degree d). If this is not the case, we proceed as follows. 

Observe that Lemma 4.16 and Proposition 4. 17(ii) imply that J2 veT> d fv{ x v) d p „ is 
in M. Subtracting this from (4.4) shows that 



^2 fv(x v )d Pv G M 

v£T\T> d 
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Thus, we can start the above procedure again with this new "smaller" vector field, each 
component being a polynomial of degree at most 

d! = max d v , 

v£T\V d 

with obviously dl < d. Defining then V d — {v 6 T | d v = d'}, we get as in 
Lemma 4.19 that (4.15) holds for all % C V d . We then proceed like this inductively, 
dealing at each step with the components of highest degree and "removing" them, until 
all the remaining components have the same degree d~ (which is equal to min.„ e 7- d v ). 
If d~ > 2 we obtain again as in Lemma 4. 19 that (4. 15) holds for all % C V d . And if 
d~ = 1, the conclusion follows from (4.12). Thus, (4.15) holds for every equivalence 
class %, regardless of the degree of the polynomials involved. The proof of Theorem 4.3 
is complete. □ 

Remark 4.20. Our method also covers the case where each particle ueV can have an 
arbitrary positive mass m v . The proofs work the same way, if we replace the functions 
/„ with /„ = V c " (x v )/(m c m v ). Thus, if for example all the V" v , v 6 T are the same, 
but all the particles in T have distinct masses, then all the new /„ are different, and the 
particles in T belong each to a separate %■ 



5. Controlling a network 

We now show how Theorem 4.3 can be used recursively to control a large class of 
networks. The idea is very simple: at each step of the recursion, we apply Theorem 4.3 
to a controllable particle (or a set of such) in order to show that some neighboring 
vertices are also controllable. Starting this procedure with the particles in V* (which are 
controllable by the definition of M), we obtain under certain conditions that the whole 
network is controllable. 

In order to make the distinction clear, we will say that a particle c is a controller if it 
is controllable and if we intend to use it as a starting point to control other particles. 

Definition 5.1. Let J be the collection of jointly controllable sets (i.e., of sets A <Z V 
such that >~2v£A ®Pv e -M, and therefore also >~2 v eA ®Qv by Lemma 3.1). 

Obviously, a particle v is controllable iff {v} e J . The next lemma shows what we 
"gain" in J when we apply Theorem 4.3 to a controller c. Remember that the set T c of 
first neighbors of c is partitioned into equivalence classes 7^ c , as discussed in Sect. 4. 

Lemma 5.2. Let c e V be a controller. Then, 

( i) for all i, 

( ii) for all i and for all A E J the sets 

An% c , A\Ti C and % C \A (5.1) 

are in J. 



We illustrate some possibilities in Fig. 5.1. 
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(a) (b) (c) 



Fig. 5.1: a: A controller c and a few sets in J, shown as ovals, b: The equivalence 
classes 7~ c are shown as rectangles. (Only the edges incident to c are shown.) c: New 
sets "appear" in J. In particular, c' and c" are controllable. 

Proof, (i). This is an immediate consequence of (4. 1) and the definition of J. 
(ii). We consider an equivalence class 7? and a set A G J. By (4.2) we find that 

[ ( qc ~ qv + 6 ov)d Pv ,^2d q J= d Pv - 1 ceA ■ d Pv 

veTf veA veAnTf veTf 

is in M.. By the linear structure of M. and since 2~2 v eT c <9 PtJ is in M. by (4.1), we can 
discard the second term and we find 2~2 v eAnT- c ®Pv e -M. This proves that A n 7~ c is 
in J. Then, subtracting J^veAnTf d Pv from HvzA d p, ( res P- from J2veT. c d Pv) shows 
that J2veA\T c dpv ( res P- J2veT c \A ®Pv) i s m which completes the proof of (ii). □ 

We can now give an algorithm that applies Lemma 5.2 recursively, and that can be 
used to show that a large variety of networks is controllable. 

Proposition 5.3. Consider the following algorithm that builds step by step a collection 
of sets W C J and a list of controllable particles K. 

Start with W = {{v} v G V*} and put (in any order) the vertices ofV* in K. 

1. Take the first unused controller c G K. 

2. Add each equivalence class to W. 

3. For each Tf and each A G W add the sets of (5.1) to W. 

4. If in 2. or 3. new singletons appear in W, add the corresponding vertices ( in any 
order) at the end of K. 

5. Consider c as used. If there is an unused controller in K, start again at 1. Else, stop. 

We have the following result: if in the end K contains all the vertices of V, then the 
network is controllable. 

Proof. By Lemma 5.2, the collection W remains at each step a subset of J, and K 
contains only controllers. Thus, the result holds by construction. □ 

The algorithm stops after at most | V| iterations, and one can show that the result does 
not depend on the order in which we use the controllers. This algorithm is probably the 
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easiest to implement, but does not give much insight into what really makes a network 
controllable in our criteria. For this reason, we now formulate a similar result in terms 
of equivalence with respect to a set of controllers, which underlines the role of the 
"cooperation" of several controllers. 

Definition 5.4. We consider a set C of controllers and denote by 7V(C) the set of first 
neighbors of C that are not themselves in C. We say that two particles v, w G A/"(C) 
are C-siblings if v and w are connected to C in exactly the same way, i.e., if for every 
c G C the edges {c, v} and {c, w} are either both present or both absent. 

Moreover, we say that v and w are C -equivalent if they are C-siblings, and if in 
addition, for each c G C that is linked to v and w, we have that v and w are equivalent 
with respect to c(i.e., there is a S G R such that V" v { ■ ) = V" w { • + S)). 

The C-equivalence classes form a refinement of the sets of C-siblings, see Fig. 5.2. 




(a) C-siblings (b) C-equivalence classes 



Fig. 5.2: Illustration of Definition 5.4. We assume that all the springs are identical, 
except for the edge {ci,v\}. The particles v\, . . . , v 6 form 4 sets of C-siblings, with 
C = {ci, C2, C3}. The one containing v\ and V2 is further split into two C-equivalence 
classes, since vi and w 2 are by assumption inequivalent with respect to c\. 



Proposition 5.5. Let C be a set of controllers. Then, for each C-equivalence class U C 
Af(C), we have U G J. 

Proof. See Fig. 5.3. Let U = {v\, . . . , v n } C M(C) be a C-equivalence class. We 
denote by c\, . . . , Cfe the controllers in C that are linked to v\, and therefore also to 
V2, ■ ■ ■ , v n , since the elements of U are C-siblings. For each j G {1, . . . , k}, there is a 
T^ 3 with vi,...,v n G if 3 , and we define Sj = T^ 3 \ C. We consider the set 

k 

J=l 

Clearly, U C U, and U G J by Lemma 5.2. We have U = {vi, . . . , v n , v\, . . . , v*}, 
where the v* are those particles that are equivalent to v\ , . . . , v n from the point of view 
of C\, . . . , Cfe, but that are also connected to some controller(s) in C \ {ci, . . . , c^}. In 
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particular, for each j e {1, . . . , r}, there is a c* e C \ {ci, . . . , Cfe} and an i such that 
v* is in 5* = 7f i . By construction, 5* n U = (j). Thus, 

r 

u\\Js* = u. 

3=1 

Starting from U £ J and removing one by one the S*, we find by Lemma 5.2(ii) that 
U is indeed in J, as we claim. □ 




Fig. 5.3: Illustration of the proof of Proposition 5.5 for identical springs. We consider 
the C-equivalence class U = {^1,^2}. where C contains ci, c 2 , c* and possibly other 
particles (not shown) that are not linked to v\, v%. With the notation of the proof, we 
have S\ = {w, Vi,V2, w*} and ^2 = {v\, V2, v*} so that U = S\ fl <S 2 = {v\,V2, v*}. 
Since vl belongs to = {to, v*}, we find U \ S* — U. 

An immediate consequence is 

Corollary 5.6. Let C be a set of controllers. If a vertex v s M(C) is alone in its C- 
equivalence class, then it is controllable. 

Applying this recursively, we obtain 
Theorem 5.7. We start with Cq = V*. For each k > 0, let 

Cfe+i = Ck U {v G M{Ck) I no of/zer verfex in M{Ck) is Ck-equivalent to v} . 
Then, if Ck = V for some k > 0, f/ze network is controllable. 

Proof. By Corollary 5.6 we have that each Ck contains only controllers (remember 
also that V* contains only controllers by the definition of M). Thus if Ck = V for some 
k > we find that all vertices are controllers, which is what we claim. □ 

6. Examples 

In this section we illustrate by several examples the possibilities and the limits of our 
controllability criteria. 

Example 6.1. A single controller c can control several particles if the interaction po- 
tentials between c and its neighbors have pairwise inequivalent second derivative. See 
Fig. 6.1. 
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Vl ^Ov 2 

y v 3 fy^-o «3 




c 



Vr. 



n 




Fig. 6.1: If no two springs are equivalent, the Vi are controllable. Springs from the Uj 
to other particles or from one Vi to another may exist but are not shown. They do not 
change the conclusion. 

The example above does not use the topology of the network (i.e., the notion of 
siblings), but only the inequivalence due to the second derivative of the potentials. We 
have the following immediate generalization, which we formulate as 

Theorem 6.2. Assume that Q is connected, that V* is not empty, and that for each v G 
V, the first neighbors of v are all pairwise inequivalent with respect to v (i.e., no two 
distinct neighbors u,w of v are such that V" u ( ■ ) = V" w ( ■ + S) for some constant 
S G R). Then, the network is controllable. 

Proof. We use Theorem 5.7. Observe that under these assumptions, we have at each 
step Ck+i = Ck U jV(Cfc). Thus, since the network is connected, there is indeed a 
k > such that C k = V. □ 

One can restate Theorem 6.2 as a genericity condition: 

Corollary 6.3. Assume that Q is connected, that V* is not empty, and that for each 
e G £ the degree of the polynomial V e is fixed (and is at least 3). Then, Q is almost 
surely controllable if we pick the coefficients of each V e at random according to a 
probability law that is absolutely continuous w.r.t. Lebesgue. 

Example 6.4. The ID chain (shown in Fig. 6.2) is controllable. Our theory applies when 
the interactions are polynomials of degree at least 3; for a somewhat different variant, 
see [5]. To apply our criteria, we start with C = {c}. Clearly, v\ is alone in its C- 
equivalence class, and is therefore controllable by Corollary 5.6. We then take C = 
{c, vi}. Since v 2 is alone in its C -equivalence class, it is also controllable. Continuing 
like this, we find that the whole chain is controllable. 



Observe that the chain described in the example above is controllable whether some 
pairs of springs are equivalent or not. There are in fact many networks that are control- 
lable thanks to their topology alone, regardless of the potentials. We have in particular 



c 




Fig. 6.2: A one-dimensional chain. 
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Example 6.5. We consider the network in Fig. 6.3(a), with the set C = {ci, . . . , C4} of 
controllers. Since no two Vi are C-siblings, they are each in a distinct C-equivalence 
class (regardless of the potentials). By Corollary 5.6 they are controllable. We use this 
argument recursively to show that the network in Fig. 6.3(b) is controllable. Starting 
with C = {ci, . . . ,04}, we obtain that the particles in the second column are also 
controllable. Then, doing the same with the set C that contains these particles and the 
Ci, we find that the third column is controllable. Continuing like this, we gain control of 
the whole network. 




Our theory is local in the sense that the central tool (Theorem 4.3) involves only 
a controller and its first neighbors. When we "walk through the graph," starting from 
V* and taking at each step control of more particles, we only look at the interaction 
potentials that involve the particles we already control and their first neighbors. We 
never look "farther" in the graph. This makes our criteria quite easy to apply, but this is 
also the main limitation of our theory, as illustrated in 

Example 6.6. We consider the network shown in Fig. 6.4, where c is a controller. If 
and V c " 3 are equivalent, then our theory fails to say anything about the controlla- 
bility of the network. In order to draw any conclusion, one has to look at "what comes 
next" in the network. Of course, if the lower branch is an exact copy of the upper one 
(i.e., if the interaction and pinning potentials are the same), then the network is truly 
uncontrollable, and this is obvious for symmetry reasons. However, without such an 
"unfortunate" symmetry, the network may still be controllable. Indeed, by the study 
above, we know that the vector field Y = d qvi + d qvs is in A4. By commuting with 
X and subtracting some contributions already in M, one easily obtains that the vector 
field 

K d P^ + K^s + Kv 2 ■ (V - d P J + K'U • (d P , 3 dp J 

is in M. Observe that now the pinning potentials U Vl and U V3 as well as the interaction 
potentials V VlV2 and V V3V4 come into play. Taking first commutators with Y and then 
taking double commutators among the obtained vector fields, one obtains further vector 
fields of the form J2i=i 9i (ivi > Qv 2 > 1v 3 , <fo 4 )d Pv . , where the gi involve derivatives and 
products of the potentials mentioned above. In many cases, these are enough to prove 
that the network in Fig. 6.4 is controllable, even though our theory fails to say so. 
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vi V V1V2 D2 



c 



Fig. 6.4: The network used in Example 6.6. If V CV3 is equivalent to V CVl our theory does 
not allow to conclude, but the network might still be controllable. 



One question that might arise is: why does only the second derivative of the interac- 
tion potentials enter the theory? The next example shows that this issue is related to the 
notion of locality mentioned above. 

Example 6.7. We consider the network in Fig. 6.5, where c is a controller. We study the 
case where 

V cv (q c - q v ) = (q c - q v ) 4 , U v (q v ) = q% , 

V cw (q c - q w ) = (q c - q w ) A + a ■ (q c - q w ), U w (q w ) = q% + b ■ q w , 

for some constants a and b. The terms in a and b act as constant forces on c and w. Since 
V" v ~ V" w , the particles v and w are equivalent with respect to cby our definition. Thus, 
our theory fails to say anything. We seem to be missing the fact that when a ^ 0, the 
particles v and w can be told apart due to the first derivative of the potentials. However, 
having a ^ is not enough; the controllability of the network also depends on b. 
Indeed, if a = b, the vector field Xq is symmetric in v and w, and therefore the network 
is genuinely uncontrollable. If now a ^ b, we have checked, by following a different 
strategy of taking commutators, that the network is controllable. Consequently, when 
two potentials have equivalent second derivative, but inequivalent first derivative, no 
conclusion can be drawn in general without knowing more about the network (here, it is 
one of the pinning potentials, but in more complex situations, it can be some subsequent 
springs). 



V 




w 



Fig. 6.5: The network discussed in Example 6.7. 



Our theory applies only to strictly anharmonic systems, since we assume that the 
interaction potentials have degree at least 3. The next example shows what can go wrong 
if we drop this assumption. Again, this is related to the locality of our criteria. 
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Example 6.8. We consider the harmonic system shown in Fig. 6.6. The vertex c is a 
controller, and all the pinning potentials are equal and harmonic, i.e., of the form Ax 2 /2. 
The interaction potentials are also harmonic. The spring {c, v{\ has coupling constant 
2, the springs {c, v 2 } and {v 2 , W3} have coupling constant 1 and the spring {v^, V4} has 
coupling k > 0. Since V" Vl = 2 and V^ 2 = 1, the particles v\ and v 2 are inequivalent 
with respect to c. Yet, this is not enough to obtain that they are controllable (unlike in the 
strictly anharmonic case covered by our theory). With standard methods for harmonic 
systems, it can be shown that the network is controllable iff k ^ 2. When k = 2, one of 
the eigenmodes decouples from the controller c, and no particle except c is controllable. 
Thus, one cannot obtain that v\ and v 2 are controllable without knowing more about 
the potentials that come farther in the graph. 

VI 

c 

VI 1 V3 k V4 

Fig. 6.6: A harmonic network that may or may not be controllable depending on the 
value of the coupling constant k. 



7. Comparison with other commutator techniques 

It is perhaps useful to compare the techniques used in this paper to those used elsewhere: 
To unify notation, we consider the hypoellipticity problem in the classical form 

L = X Q + J2xf. 

i>0 

In [5], the authors considered a chain, so that V* is just the first and the last particle in 
the chain. Starting with d Pl (the left end of the chain) one then forms (with simplified 
notation, which glosses over details which can be found in that paper) 

d qi = [d Pl , X ] , d P2 = (Mi )2 ) _1 [9 5l ,Xo] , d q2 = [d P2 , X ] , 

and so on, going through the chain. Here, the particles are allowed to move in sev- 
eral dimensions, and Mjj + \ is basically the Hessian matrix of Vjj+i. This technique 
requires that Mjj + i be invertible, which implies some restrictions on the potentials. 
Villani [12] uses another sequence of commutators: 

C = {Xi} i>0 , Cj+i = [Cj,X ] + remainder, . 

With this superficial notation, the current paper uses again a walk through the net- 
work, but the basic step involves double commutators of the form 

\Zi ■■ z 2 j 

with Zi typically of the form ^ g v (x v )d Pv , where we use abundantly that the V e are 
polynomials. This allows for the "fanning out" of Fig. 6.1 and is at the basis of our 
ability to control very general networks. In particular, this shows that networks with 
variable cross-section can be controlled. 
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A. Vandermonde determinants 

Lemma A.l. Let a, . . . ,c n e R be distinct and non-zero, and let s > 0. Then, for all 
k G {1, . . . , n} there are constants n , . . . , r n e R such that for all j — l,...,n, 



Sjk ■ 



Proof. We have that the Vandermonde determinant 



„s+l s+2 
-1 c l 
„s+l s+2 



„s+l „s+2 



• c: 



-1 



s+n 



u^ 1 



1 Cl c\ 
1 c 2 c\ ■ ■ ■ a 



,n-l 



1 r r 2 ■ ■ ■ r™ _1 



n^ +i n ^-^ 

j=i+i 



i=i 



is non-zero under our assumptions. Thus, the columns of this matrix form a basis of 
R n , which proves the lemma. □ 
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